TECHNOLOGY OVERVIEW 
BEFORE TOKEN/NETWORK: 

SUPER PRE BETA DRAFT WORK IN PROGRESS DESCRIPTION OF THE PROPRIETARY ALGORITHM. TBD MEANS: “I KNOW, BUT CAN’T TELL YOU YET.” 


DISCUSSION: 


HTTPS: //DISCORD. GG/GVYFNKU 



RICHARD HEART 

EMAIL: RICHARDHEARTAUTHOR@GMAIL.COM 


COMPUTER MODELLING MUST BRING US INTO THE 
FUTURE! 

• Most of mankind's greatest achievements were modeled in a computer before we could make them in the real world. Computer modeling is key to the 
revolution that will bring mankind into the future we've been promised in science fiction. 

• Nuclear Fusion 

• Limitless free energy from nuclear fusion will only be possible through modelling plasma fusion containment geometries. 

• Traditional power generation 

• Oil and gas, where is the best spot to drill? How deep? Reactor vessels, cooling towers, steam turbines. 

• Rockets 

• Those rockets that will be taking us to mars; they don't get built until they're modelled in the computer first. The fuel cryostorage, the reactor, the 
nozzles, the descent flaps, all get modelled in computer first. 

• Space travel 

• Even the space suit needs modelling due to the pressure differential of space's vacuum. 

• Aeronautical 

• More efficient aircraft that get passengers where they want to go, faster, and at greater fuel savings. 

• Nautical 

• More efficient hull geometries, skin texture, sail geometry, propeller geometry. The bubbles tankers blow underwater to reduce water friction, and the 
cylindrical wind turbines for power deck side. 


COMPUTER MODELLING MUST BRING US INTO THE 
FUTURE! 


• Biotech 

• Artificial hearts, joint replacements, Breathing machines, and nanobots all must be modelled for efficacy and safety before the first prototype is built. 

• Drones 

• Better rotor geometries, more aerodynamic hull structure, foul weather optimization, vastly easier in the computer than in prototyping. 

• Extinction level events 

• How can we prevent Yellowstone from bowing it's top off and causing an ice age? How will the volcano react to different interventions such as 
cooling, or gradual pressure release? 

• Global warming 

• How much of it is caused by humans, were we already headed for an ice age or a greenhouse? What interventions can have the greatest impact at 
the lowest cost. 

• Natural disasters 

• Storms have been more powerful recently than in recorded history. What affect this new activity have on rainfall patterns, and crop growth? 

• Basically, if you want to make anything in the real world that has air, water, steam, gas, aerosols, droplets or plasmas flowing around or through it, you need 
to model it in the computer first, before you spend tons of money building it. 

• Now that you know how vital computer modelling is to mankind's progress, see the problems engineers and scientists face today. 

• You've heard of super computers right? Well, super computers exist, for the most part, to perform exactly that. The giant rows of computers full of blinking 
lights are working right now, trying to find and understand solutions to all of problems you've seen above and more! 


"TBD" TOKEN SOLVES ALL THESE PROBLEMS! 


• Here's some of the problems super computers face right now: 

• 1 . Slow: Jobs take days to run. 

• 2. Jobs fail often: After you've waited days for the job to complete, you find out you 
didn't ask the question of the computer properly. 

• 3. Inaccurate software: They're not using the best software for the job. Too many shortcuts 
are taken, making estimates for things like boundary density, viscosity, etc. 

• 4. Monopoly control of the ecosystem. 1 company has 90% market share, and barely 
updates it's product for the last 1 5 years. The faster you want your simulation to run, the 
more they charge you. 


EXAMPLE PROCESSES IN MULTI-PHASE FLOWS 


• Particle/Droplet Transport 

• Minimizing surface deposition and losses 

• Gas turbine degradation 

• Ventilation and fdter design 

• Maximizing / controlling surface deposition and yeild 

• Administration of pulmonary medicines 

• Application of coatings and fdms in optical fiber mfg. 

• Particle/Droplet Dynamics (Formation and Evolution) 

• Pollution prediction and control 

• Combustion optimization and particle formation in exhausts 

• Particle manufacturing 

• Generation of Powders, Sprays, Nanophase Materials by gas-to-particle and 
droplet-to-particle processes (see next slide) 


EXAMPLE: PROCESS DESIGN ON THE PARTICLE SCALE 



Droplet to Particle Process 


EXAMPLE: PROCESS DESIGN ON THE REACTOR 
SCALE 


Starting Location Residence Time 

(x(mm),y(mm)) sec 

A) (0,9.0) 7.6 
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Calculation of Buoyant Flow in a Microparticle Reactor 
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COMPARISON OF THE "TBD” TOKEN 
MODEL TO AVAILABLE 
COMPETITORS 


/ 


KEY ELEMENTS OF ''TBD" TOKEN MODEL 




• Basic Fluid Flow (1-D, 2-D, 3-D, Viscous, Compressible) 

• Internal \ External (e.g. pipes, ducts and rooms \ airplanes, cars) 

• Incompressible \ Compressible \ Transonic \ Supersonic \ Hypersonic (Mach 0-5) 

• Laminar \ Turbulent (e.g. small to large volumes) 

• Strong Coupling to Aerosol Dynamics & Chemistry (e.g. multi -physics) 

• Multi-Phase Aerosol Behavior (1-D, 2-D, 3-D, Polydisperse) 

• Aerosol Transport 

• Aerosol Dynamics 

• Chemically Reaction (1-D, 2-D, 3-D, Multi Species) 

• Multi Species Transport (Laminar \ Turbulent) 

• Strong Coupling to Aerosol Dynamics & Chemistry 

• Fluid phase and Surface Reactions 


EXAMPLES OF COMPETITIVE TECHNOLOGY: 
AVAILABLE MULTI-PHASE MODELS IN ANSYS* 


• Discrete Phase Model (DPM) 

• Mixture Model 

• Volume of Fluid Model (VOF) 

• Eulerian Multiphase Flow Model 

• Fine Particle Model (3rd Party Add-On) 


*market leading CFD provider holding 90% of market 


AVAILABLE MULTIPHASE MODELS IN ANSYS 
DISCRETE PHASE MODEL (1 /I ) 

• Trajectories of Particles/Droplets/Bubbles are computed in a Lagrangian or 
moving reference frame 

• Particles can exchange heat, mass and momentum with the continuous gas phase. 

• Particle-Particle Interaction is neglected 

• Turbulent Dispersion can be modelled with a Stochastic Tracking or a “Particle 
Cloud” Model 

• Volume Loading: Volume Fraction < 12% 

• Particulate Loading: Low to Moderate 

• Application Examples 

• Cyclones, Spray Dryers, Particle Separation and Classification, Liquid Fuel and 
Coal Combustion 



AVAILABLE MULTIPHASE MODELS IN ANSYS 
MIXTURE MODEL (1/2) 

• Model N-Phase Flows 

• Solve Mixture Momentum equations (for mass-averaged mixture velocity) 

• Inter-phase exchange terms depend on relative (slip) velocities 

• Turbulence and Energy Equations are solved for the Mixture 

• Only one of the phase may be defined as compressible 

• Solves the transport equation of volume fraction for each secondary phase 





AVAILABLE MULTIPHASE MODELS IN ANSYS 
MIXTURE MODEL (2/2) 

• Flow Regime: Bubbly flow, droplet flow, slurry flow 

• Volume Loading: Dilute to moderate 

• Particle Loading: Low to moderate 

• Turbulence: Weak coupling between phases 

• Stokes Number: Stokes Number << 1 

• Applications: 

• Hydrocyclones, bubble column reactors, solid suspensions, gas 
sparging 



AVAILABLE MULTIPHASE MODELS IN ANSYS 
VOLUME OF FLUID (1 /2) 

• Model to track the position of the interface between two or more immiscible 
fluids. 

• A single momentum equation is solved and the resulting velocity is shared by all 
phases 

• Surface tension and wall adhesion effects can be taken into account. 

• Solves transport equation for volume fraction of each secondary phase 

• Recommended that simulation be performed in unsteady mode. 



AVAILABLE MULTIPHASE MODELS IN ANSYS 
VOLUME OF FLUID (2/2) 

• Flow Regime: Slug flow, stratified/free surface flow 

• Volume Loading: Dilute to dense 

• Particle Loading: Low to high 

• Turbulence: Weak to moderate coupling between phases 

• Stokes Number-. All ranges of Stokes number 

• Applications: 

• Large slug flows, filling, off-shore oil ta nk sloshing, boiling, 
coating. 



AVAILABLE MULTIPHASE MODELS IN ANSYS 
EULERIAN MULTIPHASE MODEL (1 /2) 

• Solves continuity, momentum and energy equations for each phase. 

• Volume fractions characterize equation set for each phase 

• Several models available to define inter-phase exchange coefficients. 

• Strong coupling makes this model more difficult to use than Mixture Model. 

• Euler Granular Option: each granular phase is treated as a distinct interpenetrating 
granular ‘fluid’. 

• Heat Transfer between n-phases: Ranz-Marshall (Euler/Euler), Gunn 
(Euler/granular) and user-defined models 


AVAILABLE MULTIPHASE MODELS IN ANSYS 
EULERIAN MULTIPHASE MODEL (2/2) 

• Applicability 

• Flow Regime: Bubbly flow, droplet flow, slurry flow, 

fluidized beds, particle-laden flows 

• Volume Loading: Dilute to dense 

• Particulate Loading: Low to high 

• Turbulence Weak to strong coupling between phases. 

• Stokes Number: All ranges of Stokes number 

• Application Examples 

• High Particle loading flows, slurry flows, sedimentation, hydro-transport, fluidized 
beds, risers, packed bed reactors 


THIRD PARTY ADD-ON 

FINE PARTICLE MODEL (FPM) OF CHIMERA TECHNOLOGIES 


• Based on a multi-modal lognormal moment model similar to that of "TBD 
token but less general. 

• Available as a user defined function (UDF) in Fluent 

• Applicability 

• Flow Regime: Fine particle aerosol flow 

• Volume Loading: Dilute 

• Particulate Loading: Low 

• Turbulence No coupling between phases. 

• Stokes Number: Zero Stokes number limit 

• Application Examples 

• Fine particle deposition, nanomaterials processing, atmospheric dispersion 



PROVEN "TBD" TOKEN APPLICATIONS 

• Prediction/Minimazation of Atmospheric Deposition 

• Prediction/Minimazation of Heat Exchanger F ouling 

• Aerosol Control in Nuclear Reactors 



• Design of Particle Characterization Equipment 

• Design of Dry Particle Inhalers 

• Design of Nebulizer Pulmonary Drug Delivery Systems 

• Design of Photocatalytic and Optical Particle Synthesis Processes 

• Design of High Surface Area Metal Catalyst Particle Synthesis Processes Design of 
Nanostructured Doped Silica Particles for Fiber Optic Preform Manufacture. 

• Design of Gas Phase Catalyst Synthesis Processes 

• Design of Carbon Nanotubes and Nano Onions Synthesis Processes 

• Design of Nanostructured Inhalable Drug and Drug Carrier Particle Synthesis 
Processes 

• Design of Nano-drug Particle Synthesis Processes 

• Design of Nano-particle Deposition Process 

• Etc. 



OVERVIEW OF SUB-MODEL VALIDATION CASES 




) 





Particle coagulation in uniform flow. Lee & Chen ,1984 (T,M), Pratsinis, 1988 (E, T), 

Lee et al,1984 (T), Freidlander, 1977 (T) 


Laminar 2-D channel flow Schlicting, 1975 (T) 

Laminar 3-D flow in a 90 degree bend Hille et al.,1985 (E, M) 

Laminar and Turbulent flow over 2-D and 3-D backsteps Srinivasan and Rubin, 2002 (E, M) 

Steady and unsteady supersonic flow in a VDC engine inlet Saunders and Keith, 1991(E, M) 

Depositional efficiency due to diffusion in channel flow Hinds, 1982 (E,T) 

Deposition due to diffusion & thermophoresis in channel flow Chang et al., 1990 (E, T) 

Heat transfer and buoyancy in a furnace reactor Brown, 1997 (M) 

Particle deposition in a cooled turbine cascade Brown et .al., 1996 (M, T) 

Condensation shock formation in wet steam nozzle Moore et al., 1973 (E, T, M) 

Argon particle nucleation in a converging nozzle Wu & Biswas, 1996 (T, M) 

Copper nanoparticle formation in a laminar flow reactor Nasibulin et al, 2001 (E) 

Methane/air combustion in 1-D flow Kee et al., 1985 (M) 

Hydrogen/air combustion kinetics in a supersonic flow Hagenmaier, 1997 (E, M) 

CsOH aerosol dynamics & transport in a laminar flow reactor Pyykonen et al., 2002 (M) 

Photocatalytic reduction of organics on a reacting surface Yamashita et al., 1997 (E) 

Unsteady turbulent flow in a dry particle inhaler Brown et al, 2003 (E) 

Thermal decomposition of Cu(acac) 2 in a laminar flow reactor Nasibulin et al, 2001 (E) 

Sulfu ric acid aerosol formation in a HSCT nozzle an d plume Miake-Lye et al., 1993 (E, T), Appleman, 1 953 ( E, T, M) 



COMPARISON OF THE "TBD" TOKEN 
METHOD TO COMPETITORS 

EFFICIENT CFD FOR SINGLE PHASE 
FLOWS: 

THE BULK FLOW SOLVER 


o 


“TYPICAL” CFD SOLVER (AS IN ANSYS) 




• Time-marching procedures requiring a transient solution to reach Steady State (time 
step limited by strict CFL conditions) 

• Boundary conditions are approximate 

• Fully explicit or only point or line implicit 

• Based on SIMPLE, SIMPLEC or other 40 year old Pressure Relaxation Techniques 

• Accelerated and Stabilized by 

• Multi-Grid 

• Matrix Preconditioners 

• Upwinding and Complex Roe / Van Leer Flux-Splitting (physically incorrect) 

• Significant Artificial Damping (physically incorrect) 

• Add-ons are numerically only loosely numerically coupled 

• Result: It takes a long time to converge to an inaccurate solution 


THE "TBD" TOKEN CFD SOLVER 

• Based on the Reduced Navier-Stokes (RNS) plane implicit spatial 
marching technique 

• Boundary conditions are physically correct 

• Requires no time-accurate transient solution to reach the steady state (no 
“CFL” time step condition) 

• Requires no artificial viscosity 

• Separate pressure and convective flux splitting 

• Heat transfer, species-chemical reactions and particle phase transport and 
dynamics are fully integrated into the solution algorithm 

• Efficient, stable and accurate 

• Result: It takes a short time to converge to an accurate solution 


]\ COMPARISON OF FLUENT/ANSYS AND "TBD" TOKEN 

FOR FLOW WITH HEAT ADDITION (CONVERGENCE AND ACCURACY) 



Code Time Required Grid Error Iterations 


SW ~ 3.1 mins 85 x 31 le-8 92 


o 


ANSYS ~ 28 mins 150 x 50 


le-6 22188 



FLOW MODEL: FULL NAVIER-STOKES (FNS) 
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REDUCING THE PROBLEM TO ALGEBRA 

• Ax = b > x = b 

75x75x75 grid (421,875 points), 6 variables 
2,531,250 unknowns, 6.4xl0 12 elements of Al 
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Direct Methods (Implicit) 


• Cramer’s Rule: 

(N+l)! 

• Gausian Elimination: 


• Singular Value Decomp. 

N 2 

• Thomas Algorithm: 

3N 

• Sparse Martix Methods: 

5N 


• Iterative Methods (Explicit) 



Operations 

N 3 Operations 

Operations 
Operations 
Operations 



LIMITATIONS FOR CFD ALGORITHMS 




• It is impossible to get a solution to a nonlinear PDE in one 
iteration 

• For a system of fluids PDEs, it is impossible to solve a whole 3D 
domain at once except for very small cases (approximately 
15x15x15 on our computers) 

• For additional equation sets, problems get worse 

•To get solutions to realistic problems have to break the system 
down 

• Equation sets 

• Multiple domains 

• Explicit schemes that don’t require matrix inversions 

• Multi Step schemes 

• Point, line or plane implicit schemes 

• In general, the more implicit you are, the more stable your 
solution method and the fewer steps to converge 
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<- Symmetry Plane 


CONVERGENCE ACCELERATING FEATURES 

• Reduced Navier-Stokes formulation with Deferred Corrector 

• Simplified Residual Space 

• Space Marching with Pressure Flux Splitting 

• Infinite time step (no CFL condition) 

• Single sweep solutions for supersonic flows. 

• Backs weeping for fast pressure relaxation 

• Equation Staggering 

• Assures diagonal dominance without preconditioning 

• Allows fewer grid points, small domains and thus reduced computation size 

• No artificial viscosity reduces stencil size and coefficeint matrix density 

• Plane Implicit 

• Local solutions are “stable” without underrelaxation 

• Solved in Delta Form 

• Multiple back-substitutions of each LU Decomposition (b=A _1 x) 

• Continuity always 100°/o satisfied 

• No Poisson Solver for Pressure Relaxation 
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TRAVELING IN RESIDUAL 
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TRAVELING IN RESIDUAL SPA 
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SIMPLIFYING RESIDUAL SP. 
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ACCURACY ENHANCING FEATURES 




• Generalized Curvelinear Coords. & GIE Source Terms 

• Arbitrary Structured Grids w/o orthogonality constraints 

• Greatly reduces grid related metric errors in solution 

• Transformed Coordinates for Transformed Velocity 

• Velocities in Conserved Variable Form (pU), Ht vs. T 

• Equation Staggering 

• 2nd order accuracy without oscillations (no odd-even decoupling) 

• Allows physically correct and consistent treatment of BCs 

• Eliminates spurious solutions due to BC influences 

• Independent Pressure and Convective Fluxes 

• Preserves mathematical characteristics of governing equations 

• Insures all boundary conditions are physically correct over all Mach numbers 

• Consistent Boundary Conditions (P and V) 

• Allows boundaries to be set close to action 

• Spacial Marching with Pressure Flux Splitting 

• Continuity always 100% satisfied 


"TBD" TOKEN MARCHING TECHNIQUE 
(MARCHING DIRECTION, l-J PLANE) 
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"TBD" TOKEN MARCHING TECHNIQUE 
(MARCHING DIRECTION, J-K PLANE) 
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CURRENT MARCHING TECHNIQUE 
(MARCHING DIRECTION, l-J PLANE) 
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NEW “VOLUME MARCHING” TECHNIQUE 
(MARCHING DIRECTION, l-J PLANE) 
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"TBD" TOKEN BOUNDARY CONDITIONS 
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WHY DEFER AXIAL DIFFUSIVITY? 


• In the Steady-State 

• "TBD" token velocity gradients (t xx t x ) << Normal velocity gradients (t 

y 

• Majority of flow domain is “Inviscid” 

• High shear only occurs locally and perpendicular to the flow direction (e.g very 
close to walls/jets) 

• Except 

• Re very small (Re < 10) 

• Normal Shocks (Momentum & Pressure forces dominate) 
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PRESSURE FLUX SPLITTING 

• Acoustic Flux vs. Convective Flux 

• Convective flux travels in flow direction 

• Acoustic flux travels along characteristics 





“DELTA” FORM 

• Ax = b > x = A _1 b 

• A(x-dx) = b - A dx 

or 

• A dx = A(x+dx)-b = lr > dx = A _1 b’ 

RHS is now the error (residual) 

• As dx > 0 (as the solution converges) 

b’ - — > 0 

• A ' 1 does not have to be exact 

• Thus a' 1 « A' 1 & dx = a -1 b’ 


ORTHOGONAL CARTESIAN 
VS. GENERALIZED COORDINATES 
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(K-DIRECTION, l-J PLANE) 



-2 -1 0 1 2 3 4 5 

1 2 3 4 5 6 




SPECIFICS OF THE "TBD" TOKEN 
DISPERSED PHASE MODEL 
FOR AEROSOLS AND AQUASOLS 
(PARTICLES, DROPLETS AND BUBBLES) 


o 



FEATURES OF "TBD" TOKEN 
DISPERSED PHASE MODEL 

• Stand Alone Code Including Multi-phase Solver and Gridding 

• Incompressible/Supersonic 

• One/Two/Three-Dimensional 

• Internal/External 

• Complex Geometries 

• Multi-Species with chemical reaction 

• Detailed particle transport (convection, diffusion, thermophoresis) 

• Detailed particle dynamics (nucleation, condensation, reaction) 

• Polydisperse aerosols from nanometer to mm size 

• Submodels for various aerosol synthesis phenomena 

• Laminar/Turbulent (Flow/Species/Aerosol) 

• Strong 3-Way Coupled Flow/Chemical Species/Particles, Bubbles 

• Accurate 

• Robust 




SOME EXAMPLES AND VALIDATION 
OF THE BULK FLOW SOLVER 
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INCOMPRESSIBLE (MACH = 0) FLOW IN A 1 80 
DEGREE TURNAROUND DUCT 



UJmm/s) 


INCOMPRESSIBLE (MACH = 0) FLOW IN A 180 
DEGREE TURNAROUND DUCT 
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FLOW AND TEMPERATURE I 
MICROPARTICLE FURNACE I 
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Comparison of Static Pressure on the VDC Centerbody 
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FLOW IN POROUS MEDIA 





Pilot Scale Ti02 Spray Calcination Reactor 
Temperature Contours Along the Center Plane 

(Flow from Bottom to Top) 

Peak Wall Temperature = 1 040 K 
Inflow Mass Flow Rate = 36 NL/min 
Inflow Velocity = 0.306 m/s 

Inlet Diameter = 5 cm 

Reactor Diameter =13 cm 

Outlet Diameter =10 cm 

Note: Radial Direction Stretched 
by a Factor of 3. 
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SCHEMATIC OF A TYPICAL INHALATION DRIVEN 
MULTIDOSE DRY POWDER INHALER 
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AEROSOLIZATION AND DEAGGLOMERATION IN A 
MODEL DRY POWDER INHALER & PHARMACOPIA 
INLET 



High Pressure Drop Inlet 
Dosing Chamber 
Exit Nozzle Throat 
Inhaler Exit 



Model I nhaler USP Inlet (horizontal section) 
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Contours of Non- Dimensional. Turbulent Viscosity in model DPVUSP wlqt at 
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HEAT EXCHANGE IN A NANOPARTICLE REACTOR 
"TBD" TOKEN CALCULATION 


Streamlines and Contours of Gas Temperature (K) 

T wall = 800 C 
Q , = 0.7 L/min 
Q co-now = 32 L/min 
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Mole Fraction of Major Species vs. Axial Location 



Mole Fraction of Minor Species vs. Axial Location 







ROCARBON PHOTOCATALYSIS 
CTOR 


ultraviolet lamp(20WX5) 

(black lamp) 


Experimental Apparatus 



20 30 40 

Axial Location (cm) 




A=5.812x10 -5 m 2 


A=4.555x1 O’ 5 m 2 


A=3.299x1 O ' 5 1 


(■♦■solution)' 


>6x10' 


OPTIMIZATION OF A NANOPARTICLE SPRAY 
COMBUSTION SYNTHESIS AND DEPOSITION 
NOZZLE 
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DEPOSITION SYSTEM 
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THE DISPERSED PHASE MODEL 
FOR MULTI-PHASE FLOW 
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/ 



TYPES OF PARTICLE/DROPLET/BUBBLE DYNAMICS 
AND TRANSPORT PROCESSES 












SCHEMATIC OF GRID SCALE OF DISPERSED PHASE 
PROCESSES OF INTEREST 

VOF, Embedded Interface Particle Tracking (Lagrangian), 


& Level Set Methods Sectional & Moment Methods 




WHEN IS AN EULERIAN FORMULATION OK? 




• When any particle having a given state (composition, liquid fraction, 
temperature . . . ) can be substituted for any other particle having the 
same state 

• When particle history is unimportant 

• e.g. it doesn’t matter if a water droplet was formed 4 billion years ago or 1 
second ago ... it has the same properties) 

• When the flow doesn’t allow multiple histories in the same physical 
location 

• Recirculation 

• Non-convective Particle Transport 

• When there are many particles 



REPRESENTATIONS OF PARTICLE SIZE DISTRIBUTIONS 
AMENABLE TO EULERIAN SOLUTIONS IN CFD 


1) Monodisperse 


4) Sectional 



3) Spline/QMOM 


5) Discrete + 1 ,2,3,4 


6) Discrete 




REPRESENTATIONS OF PARTICLE SIZE DISTRIBUTIONS 
AMENABLE TO EULERIAN SOLUTIONS IN CFD 


1) Monodisperse 


4) Sectional 



3) Spline/QMOM 


5) Discrete + 1 ,2,3,4 


6) Discrete 




MODEL MOMENTS OF AN ASSUMED SIZE 
DISTRIBUTION 

• Method 

• Assume a statistically large number of particles per volume 

• Treat each moment of an assumed particle size distribution as a transportable 
’’species” in an Eulerian reference frame 

• Advantages 

• First Moment is essentially a species conservation equation 

• No numerical diffusion 

• Small computational resources 

(6 + n s + 3 equations) 

• Disadvantages 

• Limited to the assumed distribution 

• Individual particles do not have individual histories 

• A great deal of information is retained at a dramatically reduced 
computational expense 


COMPUTATION COST VS. ACCURACY FOR EULERIAN 
LOGNORMAL & SECTIONAL MODELS 



Unimodal Lognormal (3 eqns) 

General Sectional with 24 bins (24 eqns) 
General Sectional with 3 bins (3 eqns) 



COMPUTATION COST VS. ACCURACY FOR EULERIAN 
LOGNORMAL & SECTIONAL MODELS 



■o 




° Actual Distribution 

Unimodal Lognormal (3 eqns) 

- General Sectional with 24 bins (24 eqns) 
General Sectional with 3 bins (3 eqns) 



MPTIONS OF THE MOMENT MODEL 


Aerosol/ Aquasol Volume Distribution is Lognormal. 
Particles/Droplets/Bubbles are Spherical 

• Use of basic kinetic theory 

• Use of a single parameter for particle/droplet/bubble size 







BLEM FORMULATION 


Moment form of the GDE (Polydisperse) 
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^NORMAL MOMENT MODEL 



Lognormal Moment form of GDE (i.t.o. M k ) 
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TRANSFORMED LOGNORMAL AEROSOL MOMENT 
MODEL 

• Lognormal Moment form of GDE (i. t. o. W k ) 






Condensation 
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DEFINITION OF 


CH # FOR MULTI-PHASE 
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INTEGRATED PROPERTIES OF THE DISTRIBUTION 


• Since the particle transport and dynamics equations are solved in 
lognormal moment form, integral distribution properties can be 
calculated from the total flux of moments out of the domain 


• The total flux of moment T at station / is 



Thus the integral geometric mean particle volume and standard 
deviation of the distribution exiting through as surface are 


In 2 a = — In 
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TINUOUS / DISPERSED PHASE COUPLING 


Fully coupled 

• Density, Enthalpy, Cp, Viscosity etc. 

Additive Source Terms 

• Mass / Momentum / Energy 
• Nucleation, Condensation, Chemical Reaction, T 
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WEAK COUPLING VS. STRONG COUPLING OF 
THE CONTINUOUS AND DISPERSED PHASES 


For Weak Coupling, the Bulk Fluid equations take the 
general form 
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OME VALIDATION AND EXAMPLES 
OF THE DISPERSED PHASE MOMENT 
MODEL 
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COAGULATION IN THE FREE-MOLECULAR AND 
CONTINUUM LIMITS 


Continuum Limit 
Comparison with Lee 1983 
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COAGULATION WITH PARTICLE GROWTH 
FROM FREE-MOLECULAT TO CONTINUUM LIMITS 
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SIMULTANEOUS NUCLEATION, CONDENSATION 
DYNAMICS IN A NANOPARTICLE REACTOR 


Diamftir 


Geometric Mean Aerosol 
Aerosol Number Concentration 
Standard Deviation 
Peak Saturation Ratio 
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Phase 



5E+13-) r9.0E-04 


1 -8.0E-04 


4E+1 * -7.0E-04I 


, 3E+13- - 


2E+13- 


1E+13- 


(3 eqns) 


0.141 

I tl. 1 tfl v — i- 


l.ll 

3.25 
71.9% 


t 


a 

-6.0E-04.g 


2 

4. 

-5.0E-04 <y) 


1.0E-04 q 

8 - 

> 

-3.0E-04c^ 


-8.0E-04'| 


HI.0E+00. 

1 5 


Discrete Sectional (300 eqns) 

I.0E-04 

0.147 (im 
1.2 x 10 8 #/Ncm3 
1.22 
3.26 
71.8% 


FLOW AND SPECIES CHARACTERISTICS IN 
NANOPARTICLE EMISSIONS SIMULATOR DEVICE 
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Velocity Vectors 
U =0.5m/s, U, ,=0.01 m/s 
T m =100C , T dil =0 C 

s m =i.o ,S dil =0.0 
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CONTROLLED ISOPROPYL-ALCOHOL DROPLET 
EVAPORATION IN A HEATED REACTOR FOR 
MICROPARTICLE PRODUCTION REACTOR 




Isopropyl- Alcohol Droplet Evaporation 
Liquid Phase Volume Cone. (M t ) Contours 
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Isopropyl-Alcohol Droplet Evaporation 
Gas Phase Species Mass Fraction Contours 
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Isopropyl-Alcohol Droplet Evaporation 
Geometric Mean Droplet Diameter Contours 
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of the Droplet Size Distribution 
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PARTICLE FORMATION IN A SUPERSONIC AIRCRAFT 
PLUME 
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Dynamics 
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GRID ADAPTATION IN A 
SUPERSONIC AIRCRAFT PLUME 


Grid Adaptation in the HSCT Plume and Near the Nozzle Exit 

[Every third sridpoint shown" 
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NANOPARTICLE FORMATION IN NUCLEAR REACTOR”TBD 

TOKEN VS. FLUENT WITH ADD-ON PARTICLE MODEL 

CsOH Saturated N? 
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\ CONDENSATION 



IN A WET STEAM NOZZLE 






V PARTICLE DEPOSITION IN A DIFFUSION BATTERY 





OP ARTICLE DEPOSITION IN A BEND 












\ THERMOPHORETIC AND DIFFUSIONAL DEPOSITION 
NANOPARTICLES IN A TEMPERATURE GRADIENT 



Velocity Contours (m's) 
Flow rate = 90.7 ccm 
dT/dy = 2444 K'cm 
Kt = 0.55 


Sampler inflow and outflow plane: 


Sampling channel inflow and outflow planes 


Temperature Contours (K) 
Flow rate = 90.7 ccm 
dT/dy = 2444 K'cm 
Kt = 0.55 


Sampler inflow and outflow planes- 


Sampling channel inflow and outflow planes 
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Particle Number Concentration Contours (#m 3 ) 
Flow rate = 90.7 ccm 
dT/dy = 2444 K cm 
Kt = 0.55 
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Velocity' Contours (m’s) 
Flow rate =27.2 ccm 
dT/dy = 2444 Kan 
Kt = 0.55 


S ampler inflow and outflow plane: 


Sampling channel inflow and outflow planes 
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Temperature Contours (K) 
Flow rate = 27.2 ccm 
dT/dv = 2444 K'cm 
Kt = 0.55 


S ampler inflow and outflow plane: 


Sampling channel inflow and outflow planes 


Particle Number Concentration Contours (#/m ) 
Flow rate = 27.2 ccm 
dT/dy = 2444 K/cm 
Kt = 0.55 


Sampler inflow and outflow plane: 




Sampling channel inflow and outflowplanes 
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Soot Particle Deposition in a Gas Turbine 



Without Thermophoresis 
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Contours of Normalized Droplet Number Concentration 

L =025 m. T =373 K. P =1 atm, U =0.5 m/s 
d gP =100 o g =l-0, P p =1 3 ^ 

H, Total Injection Velocity=l .0 m/s 

10000 u = 0.82 m/s 

) - 8750 v = 0.41 m/s 

wl=- 0.41 m/s 


Particle Velocity Vectore w/ Contours of Number Concentration 

L x =.025 m. T x =373 K, P x =l atm, U x =0.5 m/s 
d gp=100 M 111 - a g =1 -°= P p =1 g™ 

H Total Injection Velocitv=l. 00 nVs 

0.82 m/s 
0.41 m/s 

-0.41 m/s --k 


0.9000 

08000 

07000 

06000 

05000 

04000 

03000 


ill 








1 .41 02E-08 
1 .3161 9E-08 
1 .2221 7E-08 
1.12816E-08 
1 .0341 5E-08 
9 .40 1 33E-09 
8.4612E-09 
7.52106E-09 
6 58093E-09 
5 6408E-09 
4.70066E-09 
3-76053E-JJ9* 

*f1i8027E-09 
9 40133E-10 


METAL NANOCATALYST PRODUCTION BY 
CONDENSATION NUCLEATION 
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METAL NANOPARTICLE NUCLEATION IN A 
TRANSONIC 
DE LAVEL NOZZLE 
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Non-Dimensional Pressure Contours 
(r-axis is stretched to show detail) 
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Mach Number Contours 
(r-axis is stretched to show detail) 
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CLEATION OF NANOPARTICLES FROM 
HOT WIRE GENERATOR 


Gas Temperature 
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METAL NANO-CATALYST PARTICLE PRODUCTION BY 
THERMAL DECOMPOSITION OF METAL(ACAC) 2 



+ ethylene(CO) + C0 2 + acetone + H 2 0 + methylethylketone 
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NUCLEATION OF METAL NANOPARTICLE 
CATALYSTS FROM A HOT WIRE 
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_ Contours of Iron V apor Saturation Ratio 
; U^O.375 m/s, 1^=298 K, T wir =1500 K 
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.Contours of Iron Particle Number Cone. (#/nf) 
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; U, rf =0.375 m/s, T irf =298 K, TC =1500 K 
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GAS CONDITIONS IN A NANOMATERIAL 
SYNTHESIS REACTOR 


Contours of Gas Temperature (C) 
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Velocity Vectors Colored by Axial Velocity (id's) 


Tf^UOOC 
CO Flow Rate = 400 cm /min 
N 2 Flow Rate = 400 cm7min 
FI, Mole Fraction = 0.07 
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Contours of Iron Mole Fraction 


CO Flow Rate = 400 cm’ /min 
N 2 Flow Rate = 400 cm'Vmin 
H, Mole Fraction = 0.07 
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CATALYST NANOPARTICLE DYNAMICS 


Particle Number 
Concentration 


Particle 
Mass Fraction 


Effective Mean 
Particle Diameter 


0.05 =- 
^0.04 ^ 
| 0.03 - 

|o.02 - 
jo.oi ^ 
I 0 
cSo.01 


Contours of Particle Number Concentration (#/m ? ) 

T = 1200 C 
C0 n fiow Rate = 400 cm 3 min 
N, Flow Rate = 400 cm 3 min 
H, Mole Fraction = 0.07 

Unsaturated Walls. Inlet Mass Fraction Fe = le-25 
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CATALYTIC FORMATION AND GROWTH OF 
NANOMATERIALS 
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Contours of Particle Number Concentration (film) 

T, jmiC5 = 1200 C (CNT function A) 

CO Flow Rate = 400 cm /min 
N, Flow Rate = 400 cmVmin 
FL Mole Fraction = 0.07 

Unsaturated Walls. Inlet Mass Fraction Fe = le-25 



0.25 

Axial Location (m) 




Contours of Particle Geometric Mean Diamter (m) 

T = 1200 C (CNT function A) 

CuPlowRate = 400 cnf /min 
N, Flow Rate = 400 cm" /min 
H/ Mole Fraction = 0.07 

Unsaturated Walls, Inlet Mass Fraction Fe = le-25 
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FLOW IN A SPRAY 



Schematic of Spray Nozzle Configuration 
of Chen and Pereira (1996) 
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Fluid Velocity Vectors and Streamlines and 
Contours of Fluid Mach Number 
Isopropyl- Alcohol Spray Rig of Chen & Pereira (1996) 
Two-Dimensional Calculation (no evaporation) 


- U foo =6.65 m/s, 1^=373 K, P^l atm (Re=28662) 

Dronlets inieoted at x= 023 m D = 30 um rr =10 
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DROPLET COMBUSTION IN A SHOCK WAVE 


Conditions: 


P 1 0 atm 
1000 m/ 
1650 K 


Aerosol Properties vs. Axial Location 
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\ DROPLET COMBUSTIOL 


Mole Fraction of Major Species vs. Axial Location 
U =1,000 m/s, P =10 atm, T =1650 K 
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Mole Fraction of Minor Species vs. Axial Location 
U =1,000 m/s, P =10 atm, T =1650 K 
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PARTICLE TRANSPORT AND DEPOSITION IN A 
NANOMATERIAL SYNTHESIS REACTOR 




PARTICLE TRANSPORT AND 
NANOMATERIAL SYNTHESI 


Contours of Radial V elocity (m s) in Recirculating Region v (m s) 



Close-up of Contours of Particle Mass Fraction at Reactor Inflow 
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DEPOSITION IN A 
5 REACTOR 


- Contours of Particle Mass Fraction in Recirculating Region 


0.05 - 


0.05 


0.025 


0.0225 


0.02 


"§0.0175 


§ 0.015 

05 


0.0125 


0.01 


Close-up 


0.1 0.15 

Axial Location (m) 


Axial Location (m) 


M-(kg kg) 


7 2E-11 
6.7E-11 
6.2E-11 
5 7E-11 
5.2E-11 
4.6E-11 
4 IE-11 
3.6E-11 
3 IE-11 
2.6E-11 
2. IE-11 
1.5E-11 
1.0E-11 
5.2E-12 
0.0E+00 


(#/m 3 ) 


2.3E+0 

2.1 E+0 
1 9E+0 
1 8E+0 
1 6E+0 
1.4E+0 
1 3E+0 

1.1 E+0 
9.6E+0 
8.0E+0 
6.4E+0 
4.8E+0 
3.2E+0 
1.6E+0 
O.OE+O 




STRENGTH OF THE “TBD” TOKEN CORE ALGORITHM 


• A general RNS based CFD code for aerosol transport and dynamics 
has been developed 

• Applicable over a wide range of particle sizes from nm to 100 mm (High 
diffusion to high inertia) 

• Applicable over a wide range of Reynolds and Mach numbers 

• Applicable for strong coupling between phases 

• Applicable to arbitrary geometries in one, two and three dimensions 

• Includes a wide range of aerosol transport and dynamics mechanisms coupled 
with gas phase and surface chemistry 

• It is efficient and accurate 

• Fully validated 

• Conclusion 

• Faster, more accurate, more flexible and powerful computer modelling 
software, proven, in the real world. 

• Unless you are building advanced things, much of the above text was probably gibberish to you. If 
you didn’t understand it, you’ll probably enjoy the easier to understand token/utility/network 
whitepaper when it’s available. It will cover all the cool Cryptocurrency token things © 

• https: //discord. gg/gvyFnKu 


